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ABSTRACT 

We use a statistical sample of ~500 rich clusters taken from 72 square degrees of the Red-Sequence 
Cluster Survey (RCS-1) to study the evolution of ^30,000 red-sequence galaxies in clusters over the 
redshift range 0.35<z<0.95. We construct red-sequence luminosity functions (RSLFs) for a well- 
defined, homogeneously selected, richness limited sample. The RSLF at higher redshifts shows a 
deficit of faint red galaxies (to My > -19.7) with their numbers increasing towards the present epoch. 
This is consistent with the 'down-sizing' picture in which star-formation ended at earlier times for the 
most massive (luminous) galaxies and more recently for less massive (fainter) galaxies. We observe a 
richness dependence to the down-sizing effect in the sense that, at a given redshift, the drop-off of faint 
red galaxies is greater for poorer (less massive) clusters, suggesting that star-formation ended earlier 
for galaxies in more massive clusters. The decrease in faint red-sequence galaxies is accompanied 
by an increase in faint blue galaxies, implying that the process responsible for this evolution of faint 
galaxies is the termination of star-formation, possibly with little or no need for merging. At the bright 
end, we also see an increase in the number of blue galaxies with increasing redshift, suggesting that 
termination of star-formation in higher mass galaxies may also be an important formation mechanism 
for higher mass ellipticals. By comparing with a low-redshift Abell Cluster sample, we find that the 
down-sizing trend seen within RCS-1 has continued to the local universe. 

Subject headings: galaxies: clusters: general ~ galaxies: elliptical and lenticular, cD - galaxies: 
evolution - galaxies: luminosity function, mass function 



1. INTRODUCTION 

Clusters of galaxies are ideal laboratories for studying 
galaxy evolution since they contain many galaxies seen at 
the same epoch in close proximity. Their cores are dom- 
inated by early- type galaxies, which are the major com- 
ponent of the high mass end of the galaxy stellar mass 
function locally. There is now a good deal of evidence 
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that cluster early-type galaxies formed the bulk of their 
stars at high redshift and thereafter simply evolved pas- 
sively with little or no residual star-formation. One such 
line of evidence is the ti ght sequence they form in color- 
magn itude space (e.g., IVisvanathanI 119781 : iBower et al.l 
[1992), 'the red -sequence'. A similar red-sequence is also 
seen for early-type galaxies in the field out to at least 
z~l ([BcU ct al. 2004). Furthermore, all galaxies ap- 
pear to be divided into two distinct populations: the 
passively-evolving red-sequence and the actively star- 
forming 'blue cloud'. Only a small amount of residual 
star-formation (less than ~10% of the galaxies' past av- 
eraged star-formation rate) is necessary to move a galaxy 
from the red-sequence to the blue cloud. Therefore, 
early-type galaxies can provide unique insight into the 
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history of star-formation, as traced by objects in which 
star- formation has aheady been terminated. 

In the local universe, the probability of a galaxy be- 
longing to the red-sequence or bl ue cloud depends on its 
stellar mass and its environment (jBaldrv et al.ll2006[ ). It 
is likely that the other fundamental parameter govern- 
ing the properties of a galaxy is the epoch at which it is 
observed. Thus, in order to build a complete picture of 
galaxy evolution, we need to study the colors of galaxies 
as a function of mass (or luminosity), environment and 
redshifts. 

The classical picture for the forma tion of galaxies pro - 
poses a single 'monolithic collapse' (jEggen et al.lll96^ . 
with stars in elliptical galaxies bei ng formed in a single 
burst, thereafter evolving passively ([Partridge fc Peebles! 
[196 7: Sandagc et al. 1970J. This very simple model pre- 
dicts remarkably well many of the properties and scaling 
relations of elliptical galaxies. 

In the current hierarchical paradigm, structure forms 
in a 'bottom-up' sense, as galaxies and clusters are built 
from the merging of smaller units. Recently, there is 
growing evidence that star-formation has evolved in a 
'top-down' sense with more massive galaxies being most 
actively star-forming in the past and the bulk of the star- 
formation activity moving toward less massive galaxies 
as the universe ages. Although this seems intuitively 
at odds with hierarchical models, scenarios have been 
proposed in which sta r-formation progresses in this anti- 
hierarchical manner (|De Lucia et all 120061 ). Whereas 
previous generations of semi-analytic models in this hi- 
erarchical framework suggested that the most massive 
early -type galaxies should be younger than less mas sive 
ones (|Baugh et al.lll996l : lKauffmann fc Charlotlll998l ). in 
order to reconcile with the observed down-sizing trend, 
the prediction is now that although the most massive 
early-types assembled their mass later than lower mass 
early-types, the stellar mass has been built up through 
a series of gas-poor mergers which do not result in addi- 
tional star-formation. Hence the earlier formation times 
of the stellar populations in more massive galaxies is re- 
covered. 

Despit e numerous signs of merging in early - type galax- 
ies (e.g.. Ivan Dokkumlf2f)n5l : iTran et al.l[2005l ). it remains 
an open question how important mergers are in their 
formation and evolution. The problem of disentangling 
how a galaxy assembled its mass from how it assembled 
its stars is a difficult one. 

Several studies of field galaxies have reported this 
down -sizing or anti-hieararchical trend in star-formation 
(e.g.. 'Bell et al.' '2004'; 'Juneau et al.' '2005'; ' Faber et al.l 
l200a Bundy et al. 2005; Scarlata et al. 200'$. Initial 
results suggested that the conioving number density of 
massive early-type galaxies had evolved more than could 
be accoimted for by passive evolution alone, and that 
'dry nierging' of massive galaxies was required (|Bell et al.l 
120041: iFaber et all |2005[ ). More recently, it has been 
suggested that most, if not all, of the evolution can 
be attributed to the termination of star formation and 
pure passive evolution; and that a significant contribu- 
tion from merging is not required (jCimatti et al.l [20061 : 
IScarlata et "anl2007f ). 

In galaxy clusters, down-sizing appears to be sup- 
ported by the spectro scopic ages of red- sequence galaxies 
as a function of mass ([Nelan et al.ll2005h . In distant clus- 



ters (z~0.8), a deficit of faint red-sequence galaxies rela- 
tive to local cl usters has been claimed, in accordance with 
this picture (iDe Lucia et al.l 120041 : iTanaka et aT] l2005t 
Ide Lucia et al.ll2007D . However, all of these high-redshift 
works have found a large cluster-to-cluster scatter in 
their samples of sizes of approximately 1-10 clusters, in- 
dicating that a large, statistical sample is crucial to such 
studies. The relative contributions of passive evolution 
versus dry merging to explain the evolution of the num- 
ber density of early-type galaxies both in the field and in 
clusters is still an open question. 

In this paper we present results using the first statisti- 
cal sample of galaxy clusters drawn from a well-defined, 
wide area, homogenous survey covering a large redshift 
range, 0.35<z<0.95. We present the survey data in §2 
and detail our method for constructing composite clus- 
ters in §3. In §4 we examine the Red-Sequence Lumi- 
nosity Function and use the ratio of luminous-to-faint 
red-sequence cluster galaxies to trace its evolution with 
redshift and dependence on cluster mass. In §5 we discuss 
our results and compare with other studies of early-type 
galaxy evolution both in clusters and the field, and in 
§6 we present our conclusions. Throughout we assume a 
cosmology of Hq =70 km s"^Mpc"^ (and h = iJo/100 
km s^iMpc^i), flu =0.3 and fl/y =0.7. 

2. DATA 

The Red-Sequenc e Cluster Survey (RCS-1 
iGladders fc Ye^ [2005h is a two filter imagmg sur- 
vey covering ^^100 square degrees. It was designed to 
build a well-defined sample of galaxy clusters out to 
z^l using a highly- e f ficient color selection technique 
([Gladders fc Yeel 120001 l2005f) . This technique provides 
a photometric estin iate of the cluster re d shift, accurate 
to Az~0.05 (e.g., iBlindert et al] 120071: iGilbank et all 
[200l . A proxy for the mass of each cluster is produced 
by measuring the optical richness, which is obtained 
by calculating the amplitud e of the galaxy-cluster 
cross-correlation functio n, Bgc ([Longair fc Seldnerlll979l : 
lYee fc L6Dez-Cru3ll999[ ). We use a modified version of 
the Bgc parameter (Bgc_R, see iGladdcrs & Yce 2005), 
considering only galaxies with colors compatible with 
the red-sequence at the estimated redshift of the cluster. 
Although the uncertainty on BgcR for an individual 
cluster is relatively large (^^20-30%), the accuracy for 
ensembles is good, as demonstrated by the agreement 
between cosmological parameters derived from Bgcfl- 
selected samples of RCS-1 clusters a nd the current 
best- fit cosmology from other methods ([Gladders et al.l 
l2007f) . Extensive work is on-going to calibrate the scal- 
ing relation between Bgci? and mass from RCS clusters- 
Some e a.rly results hav e been p resented inlBlindert et al.l 
([2007f) : IGilbank et al] ([2007D : iBarrientos et all ([2007f) . 
Results from the X-ray se l ected CNOCl clus t ers ca n be 
found in[:%e_& ElliMso3 ([200l : llicirern] ([M)6h . 

We use data from the RCS-1 photometric catalogs, 
which are derived from moderate depth imaging data in 
the Rq- and z'-bands. The imaging was obtained with 
two mosaic cameras on 4-m class telescopes, CFH12K 
on the CFHT and MOSAIC-II on the CTIO Blanco 
telescope. Details of the data reduction are given in 
IGladders fc Ye^ ([20051) and we only give a brief account 
here. The survey is divided into 20 patches, each typ- 
ically around 2 degrees x 2 degrees. Object detection. 
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classification and photomet ry were performed on the im- 
ages using PPP (|Yeelll99i[ l. For each object, total mag- 
nitudes in the deeper of the two filters (usually R) were 
computed from a curve of growth analysis. Colors were 
measured using a 3" aperture or the optimal aperture 
from the curve of growth, if it is smaller. Total magni- 
tudes for the shallower filter were then calculated us- 
ing the total magnitude of the deeper filter and this 
co lor. Galactic extinct ion was corrected using the maps 
of lSchleeel et"all (floMl . The z' magnitudes are expected 
to be uncertain at the z'<0.10 level, and this is con- 
firmed by our own comparison of internal overlaps (and 
also additional photometry from other follow-up imag- 
ing) and the (R — z') c olors should be accurate to <0.03 
(|Gladders k Yeel[2005h . Throughout, we use magnitudes 
on the AB system unless otherwise specified. In this pa- 
per we only consider the 72 square deg rees selected to 
have the highest photometric quality fsee lGladders et ahl 
[20071 for details). 

3. CONSTRUCTING COMPOSITE CLUSTERS 

We constr uct compo si te clu sters following the basic 
technique of iLoh et al.l (|2007( ). We repeat the salient 
points of that work here, and add details pertinent 
to our analysis. To build our sample, we select clus- 
ters from the latest (December 2005) generation RCS- 
1 cluster catalog, with red-sequence estimated redshifts 
0.35< zrcs <0.95; richness, Bgc_R>500 (in units of {h^Q 
Mpc)^-^''); and detection significance, ctrcs >3.3-cr. A 
Bgcij=500 cluster corresponds to a velocity dispersion 
a ~400km s~^ (e.g.. IBlindert et al.ll2007f) . and our sam- 
ple extends up to Bgcfi'^2000 clusters which would cor- 
respond to a ~1200km s^^. 

For each cluster, we extract colors and magnitudes 
(binning 0.1 magnitudes in z'-band magnitude and 0.05 
magnitudes in (i? — z') color) for all galaxies within 
a radius of 0.5xr2oo, where the value of r2oo is esti- 
mated from the ric hness following the rela tion log r2oo — 
0.48 log Bgc - 0.95 (jBarkhouse et al.ll2007l . after convert- 
ing to our cosmology frorn their h = 0.5 to h = 0.7; see 
also lYee fc Ellingsonll2003[) . This corresponds to a phys- 
ical radius of between ~1.1 /i~^Mpc and 1.9/i~^Mpc for 
the range of richnesses used here. For some clusters, the 
circle defined by this radius may fall partially off the edge 
of our survey. We use detailed maps of the positions of 
all the CCD chips to calculate the fractional area lost 
due to survey geometry. We use this fraction to reject 
clusters where the fractional areal completeness is <0.6 
and to correct incomplete clusters. The total fraction of 
clusters rejected due to geometric considerations alone is 
around 20%. The rejected clusters are not a strong func- 
tion of redshift or richness, so the net effect is simply 
to reduce the total usable area of the survey by ^^20%. 
For the remaining clusters, more than 80% of the sample 
have fractional areal completenesses of >0.8. 

The data from these clusters are placed into redshift 
bins of Az=0.1. The typical uncertainty in redshift 
given the red-sequence color of a cluster varies from 
^0.04 to 0.08 with in RCS-1 fe.g.. IBlindert et al.l [20071 
iGilbank et al.l[2007h . We do not attempt to correct the 
colors of each cluster within each redshift bin to a com- 
mon redshift, since the systematic difference in color is 
smaller than the size of the random error due to the 



accuracy of measuring the position of the red-sequence 
in color-magnitude space, and doing so will in fact in- 
crease the dispersion of the composite color magnitude 
diagrams (CMDs). 

The composite CMDs so created contain cluster galax- 
ies plus contamination from background/foreground 
galaxies. We remove this contamination in a statisti- 
cal way by creating background CMDs by summing the 
data from all galaxies within RCS-1. In practice, we sum 
galaxies within each RCS-1 patch (disjoint areas of sky of 
typically ~4 deg^) and create a composite background by 
summing the background from each patch, weighted by 
the number of clusters each patch contributes to our final 
cluster sample. We do not explicitly remove clusters from 
these background fields, as the total area contributed 
by all the clusters within any patch is small (~2%) and 
masking the clusters each time is computationally expen- 
sive. We verified that performing this masking does not 
affect our final results. 

Errors are propagated through for each bin of the CMD 
accounting for the Poisson uncertainty due to cluster 
counts, background counts, and the variance from patch- 
to-patch. The number of cluster galaxies in the j'*'* bin 
of the composite cluster CMD, Ncj, is given by: 

i 

where Nij is the observed number of galaxies (cluster -|- 
field) in the j*^ bin of the i*'* cluster region {Ncj is al- 
lowed to be formally negative in our method. In practise 
we shall be binning over many bins of color and magni- 
tude in our analysis, such that the totals in these larger 
bins are always positive). The field contribution, nfj is 
given by: 

'^/j = 54rZ!^PJ^P' (2) 

where fli and flp represent the areas of the cluster re- 
gions and the p*^ patch (field) respectively. The patch 
weighting is: 



where rup is the total number of clusters in patch p 
contributing to the composite cluster. This gives more 
weight to data from those patches which contribute more 
clusters, usually due to the greater area they cover or 
greater uniformity of data. 

The error in the counts of the j*'* bin is the quadrature 
sum of cluster and field contributions 

a{N,,) = [a{n.,,r + a{n'f^r]'^\ (4) 
We denote the field uncertainty here as n'j:j since we al- 

1/2 

low for the Poisson uncertainty, nj:j , plus the patch-to- 
patch standard deviation, a^^ , in a manner akin to that 
of iLumsden et all (|1992f ) 

cr(ny^) = max{n/^^,aP^.). (5) 

The rationale for a semi-global background subtraction 
(i. e., the app r oach of Eq. 1-5) is discussed in more detail 
in ILoh et al.] (|2007f) . 
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TABLE 1 

Fitted CMR parameters for the composite clusters. 



z 


m*, 


{R- z') at m*, 


A[R - z')/dz' 


0.40 


18.940 


0.825 


-0.044 


0.50 


19.453 


0.925 


-0.046 


0.60 


19.908 


1.063 


-0.048 


0.70 


20.322 


1.243 


-0.049 


0.80 


20.675 


1.439 


-0.059 


0.90 


20.955 


1.632 


-0.063 



Note. — Columns show the average redshift of the composite 
cluster, the observed value of m* in z' (derived from the z=0.4 
fit and the passively-evolving model), the observed color at this 
magnitude and the slope of the relation. 

Once a background subtracted cluster CMD has been 
created, the next step is to fit the red-sequence. Several 
techniques were tried to fit the locus in the presence of 
the "blue cloud" galaxies. The method traditionally used 
is to fit using the biweight estimator (|Tukevl[l958D . We 
found that using this method on the whole CMD (par- 
ticularly in the highest redshift bins) caused poor fits, 
mainly due to the significant population of blue galax- 
ies at faint magnitudes which biased the fit, causing the 
relation to be much steeper than that estimated by eye. 
To circumvent this problem, we overlaid model tracks 
for the expected red-sequence and rejected galaxies bluer 
than 0.2 magnitudes in {R—z') than the expected colors. 
This effectively isolates the red-sequence (as confirmed 
by visual inspection), but still allows the slope and in- 
tercept of the relation to be fine-tuned. We also found 
that reducing the scale radius used to select red-sequence 
galaxies to 0.25xr2oo allows a cleaner fit and so we use 
this value for fitting the CMR. The results of the fit are 
given in Table [TJ 

Next, we subtract the best-fit color-magnitude rela- 
tion to leave a red-sequence which is horizontal in color- 
magnitude space, resampling the binned CMD using sine 
interpolation. The result is a red-sequence centered on 
zero color, with no magnitude dependence. We denote 
these corrected colors as [R — 2;')corr- Blue galaxies will 
have a stronger k-correction than red galaxies, and so 
we apply a differential k-correction to the blu e galaxies, 
as a function of their observed color (see e.g.. iLoh et al.l 
[200l . The primary effect of this correction is to dim the 
contribution of on-going star formation in the brightest 
blue galaxies, allowing cleaner selection of the brightest 
red-sequence galaxies. With these corrections, the result- 
ing magnitude of each galaxy is thus now more closely 
related to its stellar mass. This correction is discussed in 
more detail in ij4.3l 

Finally, to isolate red-sequence galaxies for further 
study, we choose to use only galaxies on the red side 
of the red-sequence, i.e., {R — z')cori>0. This elimi- 
nates contamination by galaxies blueward of the red- 
sequence whose magnitude errors may allow them to 
scatter onto the red-sequence. Since red-sequence galax- 
ies are the reddest normal galaxies at a given redshift, 
there should be no galaxies redward of them, after back- 
ground subtraction. We verified the CMR goodness- 
of-fit by centroiding the CMD data in color about the 
(i? — 2:')corr=0 line and applying a small sub-pixel shift, 
if necessary. Such color shifts were ^0.01. The red- 
sequence was then extracted by mirroring this distribu- 
tion about the {R — 2:')corr=0 line. To include the ef- 
fects of brightest cluster galaxies (BCGs) which are often 



found to be slightly bluer than the red-sequence, pos- 
sibly due to the effects of on-going star formation due 
to the accretion of cold gas a t the cluster center (e.g., 
iMcNamara fc 0'Connel]||l992f ) , we relax the requirement 
of (i? — 2;')corr>0 to (i?— z')corr>-0.5 for galaxics brighter 
than M* . This allows the inclusion of galaxies which are 
clearly seen to be separate from the blue galaxy popu- 
lation (after application of the differential k-correction, 
described above). 

4. RED-SEQUENCE LUMINOSITY FUNCTIONS 

We are now left with a CMD, constructed to con- 
tain only red-sequence galaxies. Red-sequence luminos- 
ity functions can be constructed by simply summing 
over the color bins. A crucial step is to understand 
the magnitude completeness of the sample. Modelling 
the incompleteness, as is usually done for galaxy number 
counts in the field, becomes a much more complicated 
problem due to the color cuts imposed. The only re- 
liable way to verify the completeness for red-sequence 
selected galaxie s is by deeper imag ing of the same ar- 
eas of sky (e.g., iCimatti et al.l 12009 ). We adopt a very 
conservative approach and cut our data to a magnitude 
limit which should provide close to 100% completeness 
for galaxies ()Yed Il99l] ). We adopt a limit 0.8 magni- 
tudes brighter than the 5-cr point source magnitude limits 
(jCladders fc Yedr2005D . This produces color and magni- 
tude limits in observed R and z' . To ensure that we are 
not incomplete for the faintest, reddest galaxies, we cal- 
culate the intercept of the i?-band magnitude limit with 
the red envelope of the red-sequence, prior to removing 
the red-sequence slope, and decrease the z' mag limit to 
this value (typically «0.2 mags brighter). This offset can 
been seen in Fig. [T] 

Since we only need the bluer (i?c-band) data to mea- 
sure the color of each galaxy, we can relax the limit from 
0.8 magnitudes brighter than the 5-cr point source limit 
{M^a,K - 0.8) to only M^„^^ - 0.3. The typical color 
errors show that we are still measuring colors of the 
faintest galaxies with an uncertainty of ^0.15 magni- 
tudes at the faintest limit adopted. We note that repeat- 
ing our analysis with the very conservative magnitude 
limits of Msct — 0.8 in both filters does not change any of 
our results at all, except to give us insufficient depth to 
make reliable measurements in our highest redshift bin 
(0.85<z<0.95), described later. 

Fig. [2] shows the red-sequence luminosity functions. 
We remind the reader that we have constructed these by 
averaging over all clusters of richness Bgci?,>500 within 
O- Sxrann. 

ICoUessI (|1989I) constructed composite luminosity func- 
tions by using: 

= — E#"' (6) 
rrij ^ Nio 

where Ncj is the number of galaxies in the jth bin of the 
composite LF, Nij is the number in the jth bin of the 
ith cluster's LF, Nio is the normalization of the ith clus- 
ter LF, rrij is the number of clusters contributing to the 
jth bin and Nco = X^i-^io- This method is optimized 
for finding the composite LF under the assumption that 
it is universal (so that a simple rescaling by the cluster 
richness is all that is needed to find the average cluster 
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Fig. 1. — Cofor-magnitude diagrams for*' the background- 
subtracted composite clusters in redshift bins of 0.1. Solid line 
indicates fit to red-sequence (Table 1). Dashed lines represent the 
100% completeness limits adopted. Note that while the complete- 
ness in z' of the RCS catalog is somewhat deeper, we adopt the 
limits marked by the dashed lines so that the catalog is complete 
in {R — z') color to colors significantly redder then the red sequence 
(see text). The 2D histograms have been Gaussian smoothed for 
display purposes. 

LF), and in the presence of cluster data extending to 
different depths for different clusters. The approach we 
adopt here is effectively setting Nio to unity, as we are 
trying to examine the average luminosity function of all 
clusters within our uniformly selected sample of clusters, 
at each redshift. By integrating over a well sampled vol- 
ume of the universe, the weighting intended by CoUess' 
Nio term actually occurs naturally, since less rich clusters 
(containing fewer galaxies) are more abundant and hence 
their contribution is 'up- weighted' relative to richer clus- 
ters. The actual number of clusters and the red-sequence 
galaxies they contribute are shown for various samples in 
Table [21 

We convert our observed z'-band magnitudes to rest 
frame V^-band, using the same method as dc Lucia ~et al.l 
(|200l . We use the GALAXEV stel lar population syn- 
thesis code (jBruzual &: Charla3l2003f ) to generate model 
galaxy SEDs arising from a single burst stellar popula- 
tion formed at =3. We use three populations of dif- 
ferent metallicities. The zeropoint of each is normalized 
such that the observed CMR of Coma (Tcrlcvich et al.l 
|2001D is reproduced. These also give reasonable agree- 
ment with our observed CMRs as a function of redshift. 
Observed magnitudes are converted to rest frame mag- 
nitudes, interpolating between the nearest m odels. We 
choos e not to explicitly fit Schechter functions ijSchechteil 
Il976l ) to each redshift bin as a) the characteristic magni- 
tude, M*, and faint end slope, a, are degenerate; and b) 
in the higher redshift bins, our data are not sufficiently 
deep to place strong constraints on the faint end slope 
using this parametric fit. Instead, to study the evolution 
of the number density of red-sequence galaxies, wc fit 
a Schechter function to the lowest redshift bin and pas- 



TABLE 2 

Number of clusters and red-sequence galaxies in each 
redshift bin for the different samples. 



BgcR>500 

main sample 



300<Bgcfl<500 
'poorer' sample 



B9cB>800 

'richer' sample 



-^clus 


^«al 


N 1 

^ " cms 




-^clus 


^«al 


0.40 57 


8378 


54 


4981 


22 


4463 


0.50 52 


7234 


75 


3942 


9 


1840 


0.60 70 


5901 


84 


3704 


16 


1682 


0.70 80 


5765 


82 


3194 


20 


2036 


0.80 103 


4253 


104 


2625 


27 


1299 


0.90 98 


3883 


63 


1165 


29 


1889 


Note. — 


N^l^s refers to 


the number 


of clusters used 



analysis, after rejecting those with too low areal completeness, etc.; 
Ngai refers to the number of red-sequence galaxies down to the 
100% completeness limit. 

sively evolve it using the above model to other redshifts. 
The data in the lowest redshift, 0.35<z<0.45, bin are 
sufficiently deep that a single Schechter function is not 
an adequate fit for the very faint est galaxies and it can 
be seen that the well-known (e.g. jBarkhouse et al.ll2007[) 
upturn occurs for dwarf galaxies at My ^-18. However, 
this fit is sufficient to be illustrative at brighter magni- 
tudes. For reference, the fitted parameters for this ref- 
erence model are a=-0.94±0.04 and m*, =18.94±0.09, 
which corresponds to an absolute M*, =-21.51. 

4.1. Luminous-to-faint ratios 

iTanaka et all (|2005l) and Ide Lucia etahl (|2007[ ) both 
used the ratio of luminous-to-faint galaxies to look for 
evolution in the faint end of the red-sequence luminos- 
ity function with redshift. Our data are not as deep 
as either of these two w orks. However , we h ave suf- 
ficient depth to use the Ide Lucia et al.l ()2007[ ) magni- 
tude limits (My > -20 and -20 < My < -18.2) to 
z=0.5, so first we compare our luminous-to-faint ratio 
with theirs. For their z=0.5 composite EDisCS cluster, 
Ide Lucia et al.l (|2007f ) find a luminous-to-faint ratio of 
0.598±0.064 (for their sample utifizing photo-z selected 
plus statistical substraction) or 0.695±0.077 (for statisti- 
cal subtraction only), measured in the observed /-band. 
We find a value of 0.644±0.029 measured in the observed 
z'-band. Our data probably more closely resemble their 
latter method, since we do not use photometric redshifts 
to reject galaxies on an individual basis. Our measure- 
ment agrees with their value to within the uncertainty. 

In order to reach a reasonable redshift to search for 
evolution in the LF, say Z'--^0.8, we must adopt a brighter 
magnitude limit. We choose to define luminous galaxies 
as -22. 7< My <-20.7 and faint as -20.7< My < -19.7. 
These limits allow us to use data up to our 100% com- 
pleteness limit at z=0.9, avoid including BCGs at the 
bright end, and provide approximately equal numbers of 
galaxies in each luminosity bin, at z=0.9. These bounds 
are shown as dashed lines in Fig. [2j The number in the 
lower right of each panel denotes the value of this ratio 
and its error. We calculate this value and its error fol- 
lowing the same method as described in §3 and this can 
be thought of as a limiting case where Ncj reduces to 
a luminosity function of two magnitude bins. We show 
the evolution of the luminous-to-faint ratio in Fig. [31 A 
simple linear fit with redshift is sufficient, oc with 
13 =0.94±0.18. 

This trend appears to be robust to changes in the exact 
choice of magnitude limits for the two bins, provided that 
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Fig. 2. — Luminosity functions for red-sequence galaxies in the composite clusters shown in Fig. [T] The mean redshift of each composite 
cluster is given in the upper left of each panel. Data are only shown down to the 100% completeness limit for galaxies. Error bars are 
dominated by uncertainty due to background subtraction (as discussed in JSj. The thick dashed line is the best-fit Schechter function of the 
lowest redshift composite cluster, evolved according to passive evolution to each redshift. The upper x-axis shows the rest-frame absolute 
V-band magnitude and vertical dashed lines indicate limits of bright and faint bins adopted in our analysis. 




Fig. 3. — Evolution of the ratio of luminous-to-faint red-sequence 
galaxies with redshift. Horizontal error bars represent redshift 
range used in each bin. Line is best linear fit accounting for er- 
rors. The open diamond shows a low-redshift comparison point 
(not include d in the fit) built from [Barkhouse et al.l 1)20071 ) data, 
discussed in i|5.1l 

the breakpoint is chosen to be My ~ —21 or fainter. 

4.2. Cluster mass dependence of luminosity function 

With our large cluster sample, we can examine the LF 
for subsamples of our data of varying richness. For this 
test we choose to split the sample into two bins using 
BgcR cuts of 300<Bgcfl<500 and Bgcfl>800. These sub- 
samples are denoted 'poorer' and 'richer' respectively. 
The typical error on BgcR is ~20%-30%, so these lim- 
its ensure that the bins are as independent as possible, 



whilst maintaining reasonable numbers in each. In Fig. [5] 
we show the evolution of the luminous-to-faint ratio, as 
in Fig. [21 but this time split by richness. The linear fit 
from Fig. [3] is overplotted as the solid line for reference. 
We see a trend in the direction that, at lower redshifts, 
richer clusters have a lower luminous-to-faint ratio than 
poorer clusters at the same redshift. At higher redshifts, 
z~0.6, the value of the ratios become indistinguishable 
within the errors, and may even reverse in the last bin. 

Incompleteness in the cluster catalog is a function of 
both redshift and richness. We examin e its effects us- 
ing completeness estimates from iGladderd (|2002i) . For 
clusters of BgcR >800, the incompleteness is negligible 
and so the richer bin is unaffected. However, for the 
poorest clusters at the highest redshift considered, the 
incompleteness can reach ~20%. To quantify the incom- 
pleteness, we integrate over the measured distribution of 
BgcR in the poorer subsample at each redshift, applying 
the completeness corrections, and compare the measured 
mean BgcR in each bin with the expected value allow- 
ing for incompleteness. To Z'^0.7 the bias in the mean 
BgcR in the poorer bin is ^10%, but by z=0.9 the clus- 
ter sample is biased 25% richer than expected. This bias 
may wash out some of the intrinsic difference between 
richer and poorer clusters at the high redshift end, if 
the difference seen at lower redshifts still exists there. 
Lowering the BgcR>300 criterion to mitigate the effect 
of this bias would likely introduce larger systematics, as 
the false positive rate is expected to significantly increase 
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Fig. 4. — As Fig. |3] but splitting the cluster sample into 
richer (Bgcij>800, filled circles) and poorer (300<Bgcii<500, open 
squares) bins. Solid line is the best fit from Fig. [3] showing the fit 
to the whole sample, and the dotted line is a linear fit to the poorer 
clusters. Dashed line shows a (l + z)*^ fit to the ric h clusters includ- 
ing a low redshift rich composite cluster based on lBarkhouse et al.l 
(|2007|) data, discussed in ^5A\ 

below this richness. 

A potentially more serious selection effect concerns the 
use of BgcR cuts. The clusters are selected by the number 
of galaxies on their red-sequence within a fixed physical 
radius. The magnitude limit adopted for the Bgcfi mea- 
surement corresponds to ^ AI*+2 or the 100% complete- 
ness limit, whichever is brighter. This is faint enough 
to be affected by the decreasing fraction of fainter red- 
sequence galaxies at higher redshifts. If clusters natu- 
rally exhibit a monotonic sequence of luminous-to-faint 
ratios which increases with decreasing richness, and we 
select clusters based on the total number of galaxies on 
the red-sequence, then this might impose a limit to the 
maximum luminous-to-faint ratio we can measure for the 
poorest clusters. This occurs since the poorest clusters 
appear to have a greater deficit of faint red members, and 
thus such systems with high luminous-to-faint ratios (low 
fractions of faint galaxies) will be systematically excluded 
from our sample. The fact that the two or three highest 
redshift bins for the poorer clusters show approximately 
constant luminous-to-faint ratios suggests that we might 
be seeing such a bias in our sample. 

Ove r the redshift range 0.4<z<0.8, Ide Lucia et al.l 
(|2007l ). when splitting their sample by velocity disper- 
sion, found a trend in the opposite direction to that 
which we see: they suggested that more massive clus- 
ters exhibited higher luminous-to-faint ratios than less 
massive clusters. We note that the 600 km s~^ division 
they used would cor respond to a richness of Bgci?.~600 
(|Blindert et al.ll2007[ ). which is very close to the dividing 
line between our richer and poorer clusters. 

4.3. Total cluster luminosity functions 

Next we consider luminosity functions for galaxies 
of all colors. In order to more fairly compare bluer 
galaxies with their red-sequence counterparts we apply 
additional corrections to the former to remove type- 
dependent star-formation differences so that the z' mag- 
nitudes more closely sample the underlying old stel- 
lar populations. Otherwise, blue galaxies temporarily 
brightened by on-going star-formation would enter our 
sample and then fade out again at lower redshift as their 
star-formation rate decreases. This is akin to deriv- 



ing a psuedo-stellar mass function, with the luminosity 
due to star formation removed. Firstly, we infer a k- 
correction by comparing the obse rved galaxy colors with 
SEDs from I Coleman et all (|1980l), shifted to the redshift 
of the composite cluster. Secondly, we add a simple 
(model-dependent but small) evolution correction to ac- 
count for the different average star-formation histories of 
the red and blue galaxies. We use a correction term of the 
form Q X z where we adopt Q=0.9 for the red-sequence 
galaxies and Q=0.6 for bluer galaxies, in the z' band. 
The Q=0.9 parameterization is an excellent approxima- 
tion to the luminosity evolution of the z=3 burst model 
we have adopted. 0=0.6 is a reaso nable choice for later 
spectral types fe.g.. lYee et al.|[2005l ). but the exact choice 
makes little (<0.1 mag) difference to the final differential 
correction over the redshift range considered here, since 
AQ<0.3 and Az=0.5. 

In Fig. [5] we plot composite luminosity functions for 
cluster galaxies of all colors. Overlaid is the same 
curve as shown in Fig. ^ showing the passively evolved 
fit to the z=0.4 red-sequence LF. This shows much 
closer agreement between the evolution expected for red- 
sequence galaxies and the observed evolution of the total 
galaxy population of all colors than with the observed 
evolution of just the red-sequence LF (c.f. data points in 
Figs. [2] and [5] with dashed line in these figures). Fig. [2] 
showed a deficit of faint red galaxies at higher redshift 
and a slow fall-off in the number of bright red galaxies 
with increasing redshift. When galaxies of all colors are 
considered, the LF (Fig. [5]) more closely resembles that 
of the passively evolved z=0.4 red-sequence LF. 

As an additional step, we examine the differences due 
to the different stellar mass-to-light ratios, M/L, of the 
red and blue galaxies by modelling the blue galaxies as 
a stellar population with an e-folding timescale of 4 Gyr. 
As pointed out by Bell et al. (2004), this color would cor- 
respond approximately to an Sb galaxy locally. This is a 
reas onable choice for a typical "blue cloud" cluster mem- 
ber (|Loh et al.ll2007D . We use the difference in M/L ratio 
between this model and that of the z=3 burst to mea- 
sure differential corrections for the bluer galaxies relative 
to the red-sequence. Applying this additional correction 
gives a stellar mass function for the total galaxies which 
still more closely resembles the expected evolution for 
the passively evolved red-sequence stellar mass function 
than does the observed evolution of the red-sequence. We 
note that our results are unchanged if instead we simply 
use the difference between the 4 Gyr e-folding model and 
the z=3 burst model to infer the differential k-, evolution, 
and M/L corrections, and use this as a check of our semi- 
empirical corrections. Given the uncertainties inherent in 
inferring star-formation histories for blue galaxies from 
a single color and modelling the entire blue cloud as a 
single, simple stellar population, we do not pursue this 
any further here, but regard these results as illustrative. 

If wc omit the corrections to the bluer galaxies just de- 
scribed and consider simply the uncorrected z'-band LF, 
then the agreement between the observed LF would even 
more closely trace the curve expected for the passively- 
evolved red-sequence shown in Fig [5] Thus, we have 
adopted a more conservative approach by applying these 
corrections instead of using the z'-band LF directly. De- 
tailed modelling of the blue populations within these 
clusters is beyond the scope of this current paper, but 
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will be explored in future work. However, this is very 
suggestive that the simplest explanation for the evolu- 
tion of the red-sequence LF is that red-sequence galaxies 
can be built up from the termination of star-formation 
in blue cloud galaxies. 

5. DISCUSSION 

We have constructed a pseudo-mass selected composite 
cluster sample from a large number of rich galaxy clusters 
out to z~l. For the first time, we have sufficient numbers 
of cluster members to study in detail the evolution of the 
luminosity function in an homogeneously selected cluster 
sample. 

5.1. Faint- end evolution 

At the faint end {My >-21), the red-sequence luminos- 
ity function in clusters declines with increasing redshift. 
This is most clearly seen in the evolution of the luminous- 
to-faint ratio (Fig. [3|) . Such a decline is consistent with 
the 'down-sizing' picture, in which star-formation pro- 
ceeds from the most massive to least massive galaxies 
as the universe ages. The red-sequence traces the his- 
tory of star-formation via the "red and dead" remnants 
of once-actively star-forming galaxies. 

In order to extend the time baseline covered by our 
RCS-1 sample, we compa re with a low red shift cluster 
sample, taking data from iBarkhouse et al.l f2007 ) who 
studied the red-sequence luminosity functions of a sam- 
ple of X-ray luminous Abell clusters in the redshift range 
0.04<z<0.20. Their method for selecting red-sequence 
galaxies was very similar to ours, except using (B — R) 
imaging (which is very close to rest-frame {R — z') at 
z^O.5). We simply convert their i?-band photometry 
into the rest frame y-band using the same stellar pop- 
ulation models described in ^ and converting their as- 
sumed h—0.5 to h=Q.7. The luminous-to-faint ratio for 
this sample is shown as the open diamond in Fig.[3l This 
low-redshift value is consistent with the extrapolation of 
the linear fit to our higher redshift data. This implies 
that, within the luminosity range we are probing, down- 
sizing is still on-going from Z'^O.Q down to the present 
day. 

5.2. Cluster selection 

In a hierarchical universe, clusters observed at some 
given epoch will have less massive progenitors at higher 
redshift. Therefore, it is pertinent to ask how the clus- 
ters at the low redshift end of our sample relate to those 
at the high red shift end. To do this, we use the results 
of Ivan den Bos ch (2002) who used N-body simulations 
to study the merger histories of cold dark matter haloes. 
From his fig. 3, it can be seen that the main progeni- 
tor of a halo of mass typical for the clusters being stud- 
ied here would grow in mass on average by a factor of 
<2 between a redshift of 0.8 and 0.4. Our mass proxy 
for selecting our sample is JigcR- Assuming the relation 
between Bgc_R and mass does not evolve between z=0.4 
and z=0.9, this would correspond to a change in the 
BgcB. limit from 500 at the low redshift end to 350 at the 
high redshift end. We have tried applying an evolving 
limit following a similar prescription and found that it 
docs not affect our main results concerni ng the evolution 
of the luminous-to-faint ratio in clusters. iGladders et al.l 



(|2007f ) found from a self-calibration technique used to 
derive cosmological parameters that the best-fit mass- 
richness relation was compatible with no redshift evo- 
lution, but with a large uncertainty. Work is on-going 
to establish whether the assumption that the Bgcj?-mass 
relation does not evolve with redshift is valid. This is 
an observationally challenging project, bu t early results 
(|Hicks et al.l 120071: iBarrientos et ail l2007t ) suggest that 
the higher redshift relation, z>0. 7-1.1, appears consis- 
tent with the lower redshift calibration adopted here. If 
anything, the relation may evolve slightly in the direc- 
tion that causes a given BgcR to represent a less mas- 
sive cluster at higher redshift, which would mean that 
the non-evolving Bgci? limit used here may naturally ac- 
count for evolution in the average cluster mass through 
hierarchical growth. 

The most serious aspect of the selection which may 
potentially influence the results is that of imposing a 
red-sequence richness cut. Bgcj? is calculated from the 
number of red-sequence members brighter than M'^-l-2 
(or the 100% completeness limit, whichever is brighter) 
within a 0.5 /i^q^ Mpc radius aperture. The coefficients 
associated with this measure are designed to make the 
value insensitive to the choice of counting radius or mag- 
nitude range under the assumption of a universal lumi- 
nosity function. However, we have shown that the faint 
end of the luminosity function for red-sequence galaxies 
both evolves relative to the bright end with redshift and 
depends on cluster richness. Below z^O.7, for the aver- 
age cluster in this study, the down-sizing effect has little 
impact on our Bgci? measurements as any effect occurs 
below the magnitude limit used for Bgcfi- Indeed, our 
red-sequence richness selection should preferentially pick 
systems least affected by down-sizing and so our results 
concerning the faint-end deficit of red-sequence galaxies 
should be regarded as a lower limit. The good agreement 
of our luminous-to-faint ratios with other work (despite 
using slightly different magnitude ranges) suggests that 
our results are not seriously affected, though. Further- 
more, the fact that mass estimates of individual clusters 
agree well wi th the expecte d ma sses based on BgcH at 
high redshift (jGilbank et al.ll2007HBarrientos et al.i2007t l 
and that cosmological constra ints based on cluster abun- 
dances ()Gladders et all l2007f ) agree with concordance 
values, suggest that this effect is likely not of serious 
concern for most of our studies. However, the cluster- 
dependent effects of down-sizing must be considered, es- 
pecially for the lowest mass systems (Fig. [4]) when study- 
ing galaxy evolution in our clusters. Accounting for these 
effects suggests a higher order correction which may im- 
prove our mass estimates based on BgcR. 

Projection of unrelated structure along the line of sight 
must also be con sidered. It is expected from simulations 
Gladder j 12002 ) and also measured from spectroscopy 
Blindert et aLlboOTtlGilbank et al.ll2007f ). that ~5-10% 



of red-sequence clusters in our study will have another 
cluster which is a significant fraction of the richness 
of the parent cluster projected onto their red-sequence. 
However, in order to be close enough in redshift space 
that the red-sequences cannot be distinguished by their 
[R—z') colors, the clusters would typically be closer than 
Az=0.1, the size of our redshift bins. So, the effect of pro- 
jections is to some extent mitigated by our bin size, as 
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the clusters are likely to be still placed in the same bin. 
The main problem is that clusters will be slightly un- 
dercounted, i.e., two clusters are placed in the same bin, 
but only counted as one. Projections where the clusters 
should be placed in different bins will just contribute a 
certain amount of cross-talk between bins, blurring the 
difference between different redshifts. However, the mag- 
nitude of this effect is likely to be ^5% from the argument 
above. Thus, the net effect of projection will mostly be 
to add noise to the RSLFs, and not to artificially pro- 
duce any of the trends we find. The richness dependence 
of the luminous-to-faint ratio means that a slight bias 
may be introduced in that the projection of two poorer 
clusters which appear as one richer cluster should have a 
slightly higher luminous-to-faint ratio than expected for 
a single richer cluster. Again, the fraction of systems so 
af fected is likely to b e small. 

ICohn et al.l |2007f ) recently used a red-sequence selec- 
tion to find clusters within the Millennium simulation. 
Their algorithm differs in detail from ours and there are 
still numerous important unresolved issues such as the 
density dependent colours of galaxies within the semi- 
analytic models used, but they estimate contamination 
rates of ~10% at z=0.4 and ~20% at z=0.9. 

5.3. Total luminosity functions 

Fig. [5] shows that, moving from high redshift towards 
the current epoch, the faint end of the red-sequence lu- 
minosity function (or, equivalcntly the stellar mass func- 
tion) becomes increasingly populated. Fig.Oshows that, 
after attempting to correct bluer galaxies' luminosities to 
stellar masses (which would otherwise allow those tem- 



porarily brightened by on-going star-formation to artifi- 
cially enter our sample), the stellar mass function of all 
galaxies in an average cluster appears relatively constant. 
Thus, the build-up of red galaxies appears to be accom- 
panied by a decline in blue galaxies in clusters. This 
suggests that a simple explanation for the build up of 
the red-sequence is the conversion of blue to red galaxies 
due to the termination of star-formation. 

The study of the evolution of the overall normalization 
of the luminosity/stellar mass function is complicated 
by the fact that our cluster mass estimator is based on 
cluster richness. Thus, evolution in number density of 
galaxies within the cluster is degenerate with evolution 
in cluster mass. However, we can use the assumption that 
the average cluster in the high redshift bins evolves into 
the average cluster in the lower redshift bins to justify 
the above argument. We showed in the previous section 
that, within the uncertainties of our sample selection, 
this is likely not a bad assumption. 

Hence, with the above caveats, we can say that not 
only do we see an increasing population of blue galaxies 
toward the faint end of the LF at higher redshift, but we 
see an increase of blue galaxies at the bright end also (c.f., 
Fig.[2]and Fig. [5]). Recent observations have found signs 
of mergers between bri ght, red galaxies in a handful of 
high redsh ift clusters (Ford et al.ll2"00l iTran et al.|[200a 
iMei et al.l [20061) and such a "dry merging" mechanism 
has been proposed for explaining the growth of the bright 
end of the red-sequence. While wc cannot place limits on 
the incidences of such mergers in our sample (due to the 
degeneracy between mass and number density) , our find- 
ing of an increase in the number of bright blue galaxies in 
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clusters at these redshifts suggests that the termination 
of star-formation may also be a significant mechanism. 
We note that although the bright end of cluster lumi- 
nosity functions locally are dominated by red-sequence 
galaxies, spectroscopically confirmed blue cluster mem- 
bers as bright as the brig htest cluster gal axies are seen 
to exist in clusters at Z'-^l (iMei et al.ll2006D. I n addition, 
in a sample of Z'^1 clustersT lFord et 311 1 20041 ) suggested 
that the brightest cluster galaxies had considerable evo- 
lution ahead of them. In 3 of their 6 clusters, the BCGs 
had morphologies of SO or later, and in at least one clus- 
ter, the BCG seemed to comprise a pair of galaxies, close 
enough to potentially merge. This is support for our find- 
ing that the red-sequence LF might be produced simply 
by conversion of star-forming galaxies to passively evolv- 
ing galaxies. 

A more thorough treatment of this problem requires 
the use of an infall model and more detailed modelling of 
the stellar populations, and such work will be presented 
in a future paper. 

5.4. Richness dependence 

Before looking at richness dependence within our own 
cluster sample, we can com pare our results with the 
COMBO- 17 field sample of ((Bell et al.|[200l . In their 
fig. 3, 1 Bell et al.l (p004;) plot a reference Schechter func- 
tion with faint end slope a =-0.6. Their actual data 
begin to fall off faster than this -0.6 line, toward higher 
redshift. In our cluster data, the faint end red-sequence 
slope never turns over as quickly as this -0.6 line, i.e., 
the field will always show a greater deficit of faint galax- 
ies in the RSLF relative to our cluster data. This sug- 
gests that star-formation ended later in faint galaxies in 
the field than in clusters. This can be viewed as an ex- 
tension of the down-sizing phenomenon: not only does 
star-formation progress from more massive to less mas- 
sive galaxies as the universe ages, but it also progresses 
from more massive to less massive clusters. 

For the cluster subsamples split by richness (Fig. H]), 
the richer sample (BgcB> 800) traces the best-fit linear 
relation (Fig. [3|) to the whole sample (Bgci?> 500), only 
with larger error bars due to fewer galaxies, again demon- 
strating that the richness cut chosen at this level does not 
affect the results. A low redshift comp arison point from 
a rich (Bgcj?>1000) subsample of the iBarkhouse et al.l 
([2007) data is also included and found to fit with the lin- 
ear relation extrapolated from our higher redshift sam- 
ple. The dashed line shows a power law fit of the form 
oc (1 + z)l^. We find a best fit /? = (1.90 ± 0.35), which 
only modestly differs from the linear fit over the red- 
shift range probed by RCS, but gives better agreement 
with the low redshift p oint from the | Barkhouse et al.l 
120071 sample. Recently, IStott et all (|2007f) parameter- 
ized the evolution in the luminous-to-faint ratio in this 
way, finding a value of /3 = (2.5 ± 0.5), which agrees 
with our value. They used an X-ray luminous sample (of 
10 clusters at z~0.1 and 10 at z~0.5 plus several addi- 
tional clusters), which should be most comparable to our 
rich subsample here. However , they probe to a fai nter 
magnitude limit (si milar to the lDe Lucia et ani2006l and 
iTanaka et al.ll2005l depths) and so it is not obvious that 
the form of the redshift evolution should be the same as 
for our sample. We note that the trend they found is in 
good qualitative agreement with our work and reasonable 



quantitative agreement. 

At the low redshift end of our RCS sample, we find that 
the poorer clusters have systematically higher luminous- 
to-faint ratios than the rich clusters, meaning that their 
faint end RSLF is falling off more rapidly than that of 
the richer clusters. This means that poorer clusters more 
closely resemble the field than do rich clusters. Again, 
none of our poorer sample's RSLF falls off as quickly 
as the up p er lim it to the field value of -0.6 found by 
iBell et al.l (|2004[ ) and thus, poorer clusters have RSLFs 
which appear intermediate between rich clusters and the 
field. Thi s is in the opposite d irection to the result 
quoted by iDe Lucia et al.l ()2006f ) in the redshift range 
0.4<z<0.8 from their sample of '^lO clusters. We sug- 
gest that the most likely cause o f this discrepancy is th e 
choice of clusters. As noted bv iDe Lucia et all (|2006f l. 
there is significant variation in the RSLFs from cluster 
to cluster, and our larger sample is more likely to be 
representative of the average cluster population. Toward 
higher redshift, we cannot measure a significant differ- 
ence between richer and poorer, given the size of our 
error bars, except in the highest redshift (z~0.9) bin 
where the trend appears to reverse. However, the last 
bin is the one which is most likely to be affected by 
incompleteness (both in the cluster sample and in the 
photometry), so we choose to be conservative and disre- 
gard this last bin. Additionally, incompleteness effects 
at high redshift act to preferentially remove poorer clus- 
ters, causing the poorer bin to shift to systematically 
higher richnesses, lowering the significance of the differ- 
ence between the two subsamples. We have overlaid a 
best-fit linear relation in Fig. |4] to guide the eye, which 
shows that a couple of points are discrepant with such a 
fit at the >l-a level. This suggests that our error bars 
may be slightly underestimated. We note that the richer 
(BgcR>800) points in Fig. [4] could be replaced with the 
main sample (Bgcfl>500) plotted in Fig. [3] to produce a 
sample with smaller error bars which would still produce 
an average composite cluster which is significantly richer 
than the poorer (300<Bgc_R<500) sample, but with more 
cross-contamination at the boundary of the richness bins, 
where the median Bgcj?, error on an individual cluster is 
ABgc_R'^200. This would still produce a significant differ- 
ence between the richer and poorer clusters at 0.4<z<0.5, 
even after increasing the size of the poorer cluster error 
bars to be consistent with the best fit straight line. 

The results presented imply that the faint end of the 
red-sequence was built up first in rich clusters, then in 
poorer clusters and finally in the field. 

6. CONCLUSIONS 

We have studied the properties of red-sequence galax- 
ies in a well-defined statistical sample of galaxy clusters 
over the redshift range 0.35<z<0.95. Each redshift bin 
of width Az=0.1 contains Ri50 clusters and «5000 red- 
sequence galaxies. Our main results are: 

1) The faint end of the red-sequence, as measured by 
the ratio of luminous-to-faint galaxies, declines with in- 
creasing redshift. This implies that star formation has 
not yet ended in the faintest cluster galaxies at high red- 
shift. The red-sequence is built up at the faint end as 
star-formation proceeds to progressively less luminous 
(less mas sive) galaxies, con sistent with the down-sizing 
scenario (iCowie et al.lfl996l) . 
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2) The turnover of the faint end of the RSLF is depen- 
dent on the cluster richness (mass) in the sense that for 
more massive clusters, the deficit of faint red-sequence 
galaxies is less than that for less massive clusters. This 
is an indication that star formation ended earlier for faint 
galaxies in richer clusters than in poorer clusters. 

3) The decline in the faint end of the red-sequence 
toward higher redshift is accompanied by an increase in 
the total (i.e., including blue galaxies) cluster LF. This 
suggests that the build up of faint, red galaxies may be 
driven largely by the termination of star-formation in 
low mass galaxies. A similar increase of blue galaxies is 
also seen at the brighter end of the LF, suggesting that 
(at least some of) the build up of high mass early-type 
galaxies may also be attributed to the termination of 
star-formation. 

Future wor k will add the B- and T^-band imaging of 
RCS-1 fields (|Hsieh et al.ll200l and the accompanying 
photometric redshift catalog to examine the luminosity 
functions of RCS clusters. The ^1000 squa re degree next 
generation survey, RCS-2 ()Yee et al.ll2007[ ). will provide 
an order of magnitude larger sample to improve upon the 
statistics of the current work. 
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